Tracking spin and charge with spectroscopy in spin-polarised ID systems 
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We calculate the spectral function of a one-dimensional strongly interacting chain of fermions, 
where the response can be well understood in terms of spinon and holon excitations. Upon in- 
creasing the spin imbalance between the spin species, we observe the single-electron response of the 
fully polarised system to emanate from the holon peak while the spinon response vanishes. For ex- 
perimental setups that probe one-dimensional properties, we propose this method as an additional 
generic tool to aid the identification of spectral structures, e.g. in ARPES measurements. We show 
that this applies even to trapped systems having cold atomic gas experiments in mind. 
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I. INTRODUCTION 
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The interest in one-dimensional structures is unwan- 
ing. For the most part, one wants to understand the spe- 
cial nature of one-dimensional materials and verify theo- 
retical models thereof. Sometimes, it is believed to help 
solving puzzles of two-dimensional systems [l| . The Lut- 
tinger liquid , as the accepted effective low-energy the- 
ory, suggests that the elementary excitations are bosonic 
collective excitations, in contrast to the quasi-particles 
of Fermi liquid theory in higher dimensions. The funda- 
mental collective excitations either carry charge ±e and 
no spin (holons/antiholons) or carry the spin ^ and no 
charge (spinous). This so-called spin-charge separation 
becomes explicit in the exact solution of the Hubbard 
model by Lieb and Wu Q using the Bcthe ansatz. 

Therefore, physicists look for experimental evidence 
of spin-charge separation in quasi-one-dimensional ma- 
terials 0. Recent promising experiments include the 
momentum-conserved tunnelling measurements between 
parallel GaAs/AlGaAs wires [6| and conductance mea- 
surements on single wall carbon nanotubes Fur- 
ther, theoretical approaches propose to use magneto- 
tunneling Q or a transport setup Q on quantum wires, 
or to extract the information from ultra-cold atomic 
gases on a lattice [loj . However, experimental techniques 
of angle resolved photoelectron spectroscopy (ARPES) 
have vastly improved over the last years, providing 
momentum-resolved spectral densities of occupied states 
in the valence bands. The most prominent experiments 
using ARPES arc those of Kim et. al. [llj on SrCu02 
and recent measurements by Claessen et. al. [l^ on the 
organic complex TTF-TCNQ, but also one-dimensional 
metal atom chains, e.g. Au on a Si(lll) substrate [l^ 
or Au on Ge(OOl) [lj| are under investigation. Every so 
often, the interpretation of experiments are critically dis- 
cussed and challenged to be ambiguous, especially when 
spin and charge scales can not be determined indepen- 
dently (see, e.g. (HI). 

ARPES experiments deliver momentum-resolved spec- 
tral densities A{k,uj) = '}2,^A'^, which probe the re- 
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being the imaginary part of the one-particle Greens func- 
tion, where H denotes the Hamiltonian of the system, 
Eq its ground state energy corresponding to the ground 

state I'I'o), are the free single particle operators cor- 
responding to a Bloch state with momentum k and spin 
a and rj is the broadening in finite systems. The spectral 
function encodes valuable information about the elemen- 
tary excitations, e.g. as a function of at fc = 0, it reveals 
the energetically separated spinon and holon excitation 
response. 

An additional homogeneous magnetic field breaks spin 
invariance and favours the population of one of the spin 
species, driving the spin degrees of freedom into align- 
ment. In this work, we propose that tuning the spin 
polarisation by such a magnetic field can be used as a 
tool to trace the spinon-like and holon-like excitations 
that are found (or claimed to be found) at zero field. We 
will show that the trajectories of the spectral features al- 
low a unique determination of the holon and thus allows 
for a definite identification of spinon and holon at zero 
magnetic field. Especially, this argument relieves experi- 
ments of the pressure to measure spin and charge energy 
scales independently. Using a density matrix renormali- 
sation group (DMRG) method we are the first to calcu- 
late the spectral function ^ in a magnetic field without 
restriction on energy or momentum. This is shown for 
the Mott-insulating Hubbard model, the Hubbard model 
in the metallic phase, and the Hubbard model in a trap- 
ping potential, comparable to cold atom gas experiments. 
As a by-product, we find uncommon signatures that shed 
light on the unexplored nature of the elementary excita- 
tions in the Hubbard model at finite magnetic. 
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II. ZERO MAGNETIC FIELD 



Starting point of our investigation is the zero-field 
spectral function of one-dimensional interacting elec- 
trons. Even for integrable models, exact analytical treat- 
ments of the spectral function are limited. For a spinful 
Luttinger liquid model, where the spin and charge sec- 
tors decouple completely, Meden and Schonhammer [Tsl 
and Voit [16[ derived algebraic singularities in A{k,ijj) 
corresponding to the holon and spinon response, which 
are, due to the nature of the Luttinger liquid limited to 
a region around uj — 0. For the simplest non-trivial, but 
experimentally relevant microscopic model, the Hubbard 
model, the celebrated Bethe ansatz offers an exact solu- 
tion by Lieb and Wu [1] and yields the collective mode's 
excitation spectrum for all values of interaction and fill- 
ing. Unfortunately, it is very tedious to extract the spec- 
tral function from the Lieb-Wu equations, but the spinon 
and holon dispersion leave a distinct footprint in the spec- 
tral function of an electron scattering state Only 
recent developments using a pseudofermion dynamical 
theory (see Ref. [l^l and references therein) facilitated 
the calculation of the spectral function in this model [l^ 
without the limitation of being perturbative in large or 
small interaction or the restriction to a low-energy effec- 
tive theory. Further, in certain Mott-insulator or charge- 
density wave insulator systems, the spectral function was 
extracted adopting a field theoretical description in the 
Luttinger framework [H, [10, EH . 

Numerical approaches are numerous, too. Exact diago- 
nalisation, as always restricted to small system sizes, was 
used e.g. on the t-J model by Eder and Ohta to com- 
pare 2D and ID spectral features. Quantum monte carlo 
methods, gen erically suffering from the sign problem, 
succeeded [2^ [13] in resolving the main spectral features 
of the Hubbard model as predicted by the Bethe ansatz 
estimates. One advantage over the DMRG method is 
the easy extension to two dimensions also studied by 
Senechal et. al. [H]. Finally, the DMRG method [H (see 
Ref. [26j for reviews) and its extensions allow for an error- 
controlled and efficient calculation of dynamical quanti- 
ties [13|- Especially Jeckelmann and co-workers did ex- 
tensive and accurate studies of the spectral properties of 
half- filled Hubbard model [1^, the Hubbard model off 
half-filling [2§|, and the extended Hubbard model (with 
nearest neighbour interaction) [sO] • The spectral features 
of TTF-TCNQ [H], e.g., seem to fit on the simple one 
band Hubbard model, which was numerically established 
by Benthien et. al. [29| and lately also by Bozi et. al. 
[18| using the pseudofermion method, but fail completely 
with a band description using density functional theory 
[l^ . Relevant to the upcoming discussion is that around 
the F-point the spectral response A{uj) displays (beside 
an incoherent continuum of spinon-holon excitations) a 
two-peak structure corresponding to holons at lower and 
spinous at higher energies. 



III. FINITE MAGNETIC FIELD 

The Hubbard model in a magnetic field was analysed 
as early as 1990 by Frahm and Korepin [3l|, [l^l with the 
focus on the asymptotics of the Greens functions in time 
and space. Even if we know of no exhaustive treatment 
for A{k,uj) for the polarised Hubbard model or other 
integrable models (see [11]), at least the low-energy be- 
havior for all momenta (corresponding to a momentum 
distribution curve in ARPES) was considered using a 
pseudofermion formulation [34| and Ref. ^33] estimated 
the behavior in the large U limit for special values of 
Lj and k. Although a recent low-energy effective field 
theory by Zhao and Liu (ssj with attractive interaction 
and spin imbalance aims at a different direction (Fulde- 
Ferrell-Larkin-Ovchinnikov state vs. Fermi liquid), they 
note that the polarisation introduces terms that couple 
the spin and the charge sector. Finally, the idea of de- 
tecting spin-charge separation by splitting the spectral 
response with a magnetic field was also formulated within 
Luttinger liquid theory by Rabello and Si [H] . Being re- 
stricted to low- lying excitations, their proposal requires 
an extreme resolution to resolve the splitting of spectral 
features. 

In contrast, we can calculate the spectral response 
A{k,u}) at any energy, momentum and interaction 
strength. As proposed, adding an external homogeneous 
magnetic field will polarise the system. At zero tem- 
perature there will be a critical magnetic field i?cr, above 
which the minority spin species will not be occupied. As- 
suming this regime is accessible in experiment, instead of 
the two peaks, one will recover the quasi-particle ((5-like) 
response, since the majority spin type now represents a 
non-interacting electron gas. Thus, we can freely tune 
between a free and an interacting electron gas response. 
Going from full to zero polarisation by turning down the 
magnetic field, we do not a priori know how the spectral 
response of an electron relates to the spectral responses of 
the charge and spin sector. Is there a smooth transition? 
Does the electronic peak split up into charge and spin? 
Or does it vanish, while the spinon and holon features 
turn up completely unrelated? 

IV. SPIN POLARISATION IN THE HUBBARD 
MODEL 

To answer this question, we calculate the spectral func- 
tion ^ for the Hubbard model 

H = -t ^ (4 ,,c^_i_^ + h.c.) -FC/^n^jn^4, 

where U = At is the Hubbard interaction strength, Cx}a 
are the local single particle fermion operators and rix.a 
is the local density of spin type a. We use a lattice 
with hard-wall boundary conditions on M sites with 
N ~ N-\ + N I electrons. For the evaluation of the resol- 
vent we employ a preconditioned Krylov-based correction 
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vector method [13, H^- We use the single particle op- 
erators corresponding to the particle-in-a-box solutions 

'^k^'^ ^ \J W+T sv[i{kx)c^,a, with k = 7:1/ {AI + G 
{!,... M}, which are known to work well for finite sys- 
tems [2^ and are equivalent to the Bloch state operators 
in the thermodynamic limit. Wc will calculate the energy 
distribution curve ^(w) for the smallest value possible in 
this expansion k = Tr/{M -1-1), since there the energy 
scales of spinon and holon are maximally spearated. En- 
ergy is measured in units of the hopping parameter t. The 
artificial broadening 77 in ([1]) is chosen to be 0.1, which 
is larger than the level spacing but small enough to re- 
solve holon and spinon features. The resolution in the 
energy is Aw = 0.05. One can attain the limit ri ^ +0 
either by a numerically unstable deconvolution [39| or by 
treating 77 as a part of the self-energy and subtracting 
it from the imaginary part of the inverse of the Greens 
function Here we use the latter. For the figures, a b- 
splined curve was used to resample to a higher resolution 
in the frequency. The DMRG cut-off is determined by a 
desired constant discarded entropy during truncation of 
0.001 and up to 1400 density eigenstates are used. 



1.2 




polarization p 

FIG. 1: Dependence of Bz in units of the gyromagnetic ratio 
7 on the polarisation in the U = 4t Hubbard model for density 
n = 26/32. 



The spin polarisation is given hy p ~ jyj -i-Af| ' ^^'^ 
instead of adding an external magnetic field B to we 
calculate in a canonical ensemble and change the number 
of up (N^) and down electrons iVj^. From the ground state 
energy for each polarisation we can derive the magnetic 
field energy contribution jSzBz = i?Gs(p) — Eqs{p = 0). 
Using Sz = \pN , we show the dependence of Bz on 
the polarisation in Fig. [1] as a step function indicating 
the range of possible magnetic field strengths for each 
polarisation. 

Examining the metallic phase of the Hubbard model 
for M — 32 sites, wc keep the density fixed at n = n| 4- 
ni = 26/32 = 0.8125, while the polarisation is stepwise 
tuned from to 1 by increasing iV| and decreasing Ni. 
Fig. [2] shows an intensity plot of 'KA{k « 0,w), where 
the ry-broadcned Lorcntz peak was deconvoluted in a way 
discussed above, retaining a residual width of 0.02 for 
visibility. 
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FIG. 2: (Colour online) Single particle spectral function 
■KA'^{k ~ 0,0;) of the U = At Hubbard model for different 
polarisations as a colour-coded intensity plot, here for the 
metallic system with average density n = 26/32. 

V. RESULTS 

On the one hand, for full polarisation we recover the 
(5-like single particle response centered at the excitation 
energy e(k « 0) = —2t. Upon turning down the magnetic 
field, the single particle peak broadens while decreasing 
in height and shifting in energy. For zero polarisation, 
we recover the established holon peak. Thus, the trajec- 
tory for the holon-like excitation is smoothly transform- 
ing from the single electron peak into the holon peak. 
On the other hand, the spinon peak at p = has a more 
interesting behavior upon increasing the polarisation. It 
has a higher weight than the holon peak at p = 0, but it 
shifts in energy and loses weight in accordance with the 
holon-like peak gaining weight. While the shoulder left of 
the holon-like trajectory can be attributed to the shadow 
band response due to k not strictly zero, there is an un- 
expected bifurcation in the spinon-like trajectory giving 
rise to a second spinon-like excitation at lower energies 
and finite field. Its weight vanishes at some finite polar- 
isation while the main excitation's weight vanishes only 
above the critical field. Also, the spinon-holon contin- 
uum at intermediate energies develops a shoulder away 
from the holon-like trajectory, vanishing together with 
and at the spinon-like trajectory. Note that we find the 
same general behavior for the half-filled Hubbard model, 
a Mott-insulator (Fig. [3]). We checked for a particular 
polarisation for M = 32, 64, 96 sized systems that these 
general features are not finite size effects. 

VI. HUBBARD MODEL IN A TRAP 

Wc add a harmonic trapping potential with Vpot > 
centered at xq = M/2 to the Hubbard model, having 
then 

The motivation is rooted in the anticipated modelling of 
a Hubbard lattice model with ultra-cold fermion gases, 
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FIG. 3: (Colour online) Single particle spectral function as 
in Fig. [5] but for the Mott-insulating state with density n = 1. 




energy 



FIG. 4: Single particle spectral function in the [/ = 4 Hub- 
bard model in a trap for strength Vpot ~ 0.01 for density 
n = 11/50 for different polarisations as a profile plot. 



which might involve such a trap. The single particle wave 
functions of the non-interacting problem (U = 0) for the 
energetically lowest states are well described by Hermite 
functions, if the trap is not too deep (V"pot(^^/2)^ ^ t). 
This leads us to use the operators 



as single particle operators in the spectral fimction ([T]), 
where g is a potential dependent normalisation constant 
and Hm {x) is the mth Hermite polynomial. The level m 
now takes the role of the momentum k in yl(fc,a;). We 
used a system of M — 50 sites with N = 22, Vpot — 0.01 
and tested for the removal of the lowest single particle 
eigenstate to = with otherwise identical parameters. In 
Fig. |4]we have plotted the spectral function 7r^(0, w) for 



selected polarisations. Even if the terminology "spinon" 
and "holons" does not apply in the trapped system, we 
recover the same overall picture of emerging features as 
in the plain Hubbard model. The bifurcation exists as 
well, but can be seen more clearly when looking at the 
polarised response (not shown) . Finally, preliminary 
calculations on an extended Hubbard model with param- 
eters as in [s^l show again a similar overall picture. 



VII. CONCLUSION AND APPLICATION 

The main message of our result is that a magnetic 
field can be used to uniquely distinguish the holon and 
the spinon peak of the Hubbard model and thus iden- 
tify the spin-charge separation. Measuring spin and 
charge energy scales independently is not necessary in 
this szenario. From zero to full spin polarisation, the 
spectral features of the charge sector continuously trans- 
forms into the electronic excitation spectrum, while spec- 
tral features of the spin sectors stay separate and vanish. 
There are, however, more facets. Firstly, our results for 
the trapped Hubbard model can encourage experimental- 
ist to realise the Hubbard model on a lattice with ultra- 
cold gases, since even with a trap, the basic notions of 
spinon and holon physics remain. Secondly, in view of 
the results for the trapped Hubbard model and the ex- 
tended Hubbard model, we conjecture that this behavior 
is generic for a spinful Luttinger Liquid. For integrable 
models this may be even analytically accessible, at least 
by looking at the excitation spectrum for finite magnetic 
field. Surprisingly, we found a bifurcation in the spinon- 
like response at a finite magnetic field, suggesting an un- 
derlying elementary excitation structure that is different 
from a simple spinon- holon picture. Further numerical 
results will be published elsewhere while, again, Bethe 
ansatz or other analytical tools may be able to reveal the 
nature of our finding. Finally, in view of experiments on 
2D high-temperature superconductor materials [l| try- 
ing to extract ID interacting properties and in view of 
the rising interest in the FFLO state [s^, [4l| for attrac- 
tive interaction and imbalanced systems, we find a broad 
range of applications and areas that may benefit from 
our proposal. 
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